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Abstract. 

In studies of cosmic -ray (CR) propagation and diffuse con- 
tinuum 7-ray emission from the Galaxy it has usually been 
assumed that the source function can be taken as smooth 
and time-independent. However, especially for electrons at 
high energies where energy losses are rapid, the effect of 
the stochastic nature of the sources becomes apparent and 
indeed has been invoked to explain the GeV excess in the 
diffuse emission observed by EGRET. In order to address 
this problem in detail a model with explicit time-dependence 
and a stochastic SNR population has been developed, which 
follows the propagation in three dimensions. The results in- 
dicate that although the inhomogeneities are large they are 
insufficient to easily explain the GeV excess. However the 
fluctuations should show up in the gamma-ray distribution 
at high energies and this should be observable with GLAST. 
Estimates of the TeV continuum emission from the plane are 
consistent with the Whipple upper limit. 



1 Introduction 

The diffuse continuum 7-ray emission from the Galaxy mea- 
sured by EGRET and COMPTEL has been the subject of 
many studies relating to CR origin and propagation. Usually 
it has been assumed that the source function can be taken 
as smooth and time-independent, an approximation justified 
by the long residence time of CR in the Galaxy, but which 
ignores the stochastic nature of the sources. 

One motivation for studying the high-energy electrons is 
the observation of the > 1 GeV excess in the EGRET spec- 
trum of the Galactic emission, which has been proposed to 
originate in inverse-Compton emission from a hard electron 
spectrum; this hypothesis can only be reconciled with the 
local directly-observed steep electron spectrum if there are 
large spatial variations which make the spectrum in our local 
region unrepresentative of the large-scale average. Strong et 
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al. (2000) presented a study of diffuse 7-rays based on a 
smooth time-independent 2D model. The new 3D model and 
its application to nuclei is presented in Strong & Moskalenko 
"New developments in the GALPROP CR propagation 
model", these proceedings (hereineafter Paper I). Here we 
concentrate on electrons and diffuse Galactic 7-rays . 



2 Effect of SNR on the electron spectrum 

The current version of the model is described in Paper I. The 
main enhancement is the inclusion of stochastic SNR events 
as sources of cosmic rays. The SNR are characterized by the 
mean time tsNR between events in a 1 kpc^ unit volume, 
and the time tcR during which an SNR actively produces 
CR. For high-energy electrons (TeV) which lose energy on 
timescales of 10^ years the effect is a very inhomogeneous 
distribution. The amplitude of the fluctuations depends on 
these two parameters, both of which are poorly known. tsNR 
is adjusted to be consistent with estimates of the SNR rate 
(e.g., Dragicevich et al. 1999). 

We consider a model with a hard electron injection spec- 
trum (index 1.8) as used by Strong et al. (2000) following 
the suggestion of Pohl and Esposito (1998) that fluctuations 
could allow such a spectrum to be consistent both with the 
GeV 7-ray excess and the locally observed electron spec- 
trum. 

Figs. 1 and 2 show the distribution of electrons with 1, 16, 
130 and 1000 GeV, for tsNR = 10'* years (corresponding to 
a "standard" Galactic SN rate 3/century). At GeV energies 
the distribution shows only small fluctuations, the particle 
density being dominated by the long storage times. At higher 
energies the losses increase and the fluctuations become sig- 
nificant as the individual SNR events leave their imprint on 
the distribution. The TeV electron distribution is quite in- 
homogeneous, but still none of the spectra around R — Rq 
resemble even remotely that observed locally (Fig. 2). For 
tsNR — 10^ years (Galactic SN rate 0. 3/century) the distri- 
bution above 100 GeV is even more inhomogeneous and the 
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Fig. 1. Distribution of 1, 16, and 130 GeV electrons at ^ = 0. 



spectrum fluctuates even more (Fig. 3). Some of the spectra 
resemble that observed locally within a factor of a few, al- 
though still none is fully compatible with the local spectrum. 

We conclude that the "hard electron spectrum" hypothe- 
sis for the EGRET 7-ray excess would require an SN rate 
much lower than standard, with correspondingly large power 
requirements for acceleration of electrons per SNR. It is pos- 
sible that the rate of CR-producing SNR could be lower than 
that of all SNR, so that a sufficiently low rate could be pos- 
sible, but this is unlikely in view of the power requirements. 
The CR luminosity (electrons + nuclei) in this model is 2 x 
lO^'i erg s^i; for 3 SNR/century we require 2 x 10^° erg/SNR 
in CR which is plausible (for 10^^ erg and 20% acceleration 
efficiency: see Fields et al. 2001 for discussion) but for 0.3 
SNR/century we need 2 x 10^^ erg/SNR in CR which would 
require quite another class of objects.' It would also result 



'We note that hypemovae may have energies of > 10^^ erg 
(Nomoto et al., 2001). 



in large fluctuations in the 7-ray emission which are not 
evident in the EGRET data (although this aspect should be 
investigated in detail). 

Our conclusion differs from that of Pohl and Esposito 
(1998), who found that a hard electron spectrum model is 
consistent with observations considering the fluctuations, but 
they included however a dispersion in the electron injection 
spectral index which increases the variations further; how- 
ever our tests with such a dispersion do not suggest that it 
eases the problem much since a large number of SNR still 
contribute to the local electron spectrum at 1 TeV. 

2. 1 Effect on the 7-ray sky 

The inverse-Compton emission becomes increasingly clumpy 
at high energies due to the effect of SNR, even for the stan- 
dard SNR rate, as shown in longitude distributions obtained 
from our model (Fig. 4). The effect is already visible at 
1 GeV and will be an important signature for the GLAST 
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Fig. 2. Distribution of 1 TeV electrons at 2 = (top) and spectral 
variations in 4 < i? < 10 kpc (bottom) for tsNR = 10^ yr. Data 
points: locally measured electron spectra; for references see Strong 
et al. (2000), with additional data from Kobayashi et al. (1999). 

7-ray observatory, which will measure up to 100 GeV. 

3 Galactic diffuse TeV 7-rays 

Recently observations of the Galactic plane (around / = 40°) 
have been reported by the Whipple Observatory (LeBohec 
etal., 2000), which place limits on the > 500 GeV intensity. 
We have extended our predicted spectrum for a hard electron 
injection spectrum to the TeV range (Fig. 5). Since the maxi- 
mum energy of accelerated electrons is unknown we consider 
the extreme case of 100 TeV. 

Even for this case the predicted spectrum from inverse- 
Compton emission is compatible with the Whipple upper 
limit, and a lower cutoff energy will be also consistent with 
Whipple; it is clear that an improved limit would quickly 
provide a critical test for models with hard electron injection 
spectra. The SNR shock-acceleration models of Baring et al. 
(1999) suggest a cutoff around 1 TeV which imply a cutoff 
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Fig. 3. Distribution of 1 TeV electrons at 2: = (top) and spectral 
variations in 4 < < 10 Icpc (bottom) for tsNR ~ 10'' yr. Data 
as Fig. 2. 

in the 7-rays around 10 GeV, in which case the predicted 
intensities are well below the Whipple limit and detection of 
TeV diffuse emission will be difficult. 
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Fig. 4. Inverse Compton 7-ray longitude distributions for a model 

with tsNR = 10* yr and hard electron injection spectrum (index 
1.8) at 6 = for 7-ray energies from 1 MeV (bottom) to 1 TeV 
(top). 



Fig. 5. Spectiiim of the inner Galaxy (330° < I < 30°, |6| < 
5°) for a model with tsNR = lO** yr and hard electron injection 
spectrum (index 1.8) and cutoff at 100 TeV. Data: EGRET (Strong 
and Mattox, 1996), COMPTEL (Strong et al., 1999), OSSE (Kinzer 
et al., 1999), Whipple (LeBohec etal., 2000). Note tiiat the Whipple 
upper limit (plotted at > 500 GeV) is for Z = 40°. 
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